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Abstract 

We propose an extension of the Landauer-Buttiker scattering theory to include effects of inter- 
action in the active region of a mesoscopic conductor structure. The current expression obtained 
coincides with those derived by different methods. A new general expression for the noise is also es- 
tablished. These expressions are then discussed in the case of strongly sequential tunneling through 
a double-barrier resonant tunneling structure. 

PACS numbers: 73.23.-b, 73.40. Gk, 73.50.Bk, 73.50.Td 



I. INTRODUCTION 



Study of electronic transport through a mesoscopic system was pioneered in the full 
quantum mechanical framework by Landauer for a single conduction channel,- and then de- 
veloped by Buttiker for multi-channel and multi-probe systems using the scattering matrix 
approach.- 1 ^ The theory, now known as the Landauer-Buttiker scattering theory (or scatter- 
ing theory, for short), has become one of the standard ways to treat the problem of transport 
through nano-scale electronic devices, which consist of at least two terminals connected to 
a domain called scattering or active region. The heart of this theory is the assumption that 
the scattering region only separates the incident wave into transmitted and reflected parts. 
Consequently, the structure is characterized by a unique set of transmission {T n } (reflection 
{Rn}) coefficients which satisfy the basic probability conservation law, T n + R n = 1. The 
transport characteristics of the structure are then expressed in terms of these coefficients 
(see Ref.jfl] for an overview). The scattering theory has successfully described the steady- 
state transport through mesoscopic structures in the ballistic regime. However, in most 
electronic semiconductor devices scattering of carriers on impurities or/and lattice vibra- 
tions is usually important, and may even govern the transport process. To solve scattering 
problems different approaches, such as the non-equilibrium Green's function^^^ 1 ^ or/and 
two-particle Green's function,—^ have been developed. Unfortunately, such approaches are 
often complicated and computationally expensive that it is not easy to gain a quantitative 
description of interaction effects. 

In contrast, the approach based on the scattering matrix is very simple in computation 
and transparent in physics. Recently, it has been efficiently developed to include the effect 
of incoherent transport on current- voltage characteristics in nano devices.— To capture such 
effects, the scattering method has been extended by considering each scattering center as 
a fictitious probe. The probe is then characterized by a phenomenological parameter in- 
terpreted as the local chemical potential, which determines a quasi-equilibrium state of the 
probe. The conservation of current at each fictitious probe, which is used to determine the 
local chemical potential, implies that when an electron escapes into the probe, another elec- 
tron has to be reinjected immediately therefrom to the conductor in such a way that there 
is no phase correlation between two events. 14 Actually, the scattering region, which includes 
a lot of scattering centers, can play a role like an especial intermediate source or sink of 
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particles. Similar to the real sources/sinks (the contacts), the current in the intermediate 
region is conserved globally rather than locally. Moreover, there should be a phase correla- 
tion between the incoming and outgoing waves. Therefore, a suitable description is required 
to feature such an intermediate region. In this work, we develop the scattering method in 
such a way that it is able to distinguish appropriately the fluxes of coherent and incoherent 
particles passing through the scattering region. The approach is then applied to discuss the 
current- volt age (I — V) and current noise (S — V) characteristics in double-barrier resonant 
tunneling structures (DBRTS). It results in a current formula coinciding with that available 
in literature, and a new general formula for the current noise power. From these formulae 
we can recover all available conclusions of the current and noise from the limit of purely 
coherent to purely sequential tunneling regime. 

The paper is organized as follows: Section II presents our key point for the development 
of the scattering theory. A current formula is then derived. In Sec. Ill, we derive a general 
formula for the shot noise power. We discuss in Sec. IV some important aspects resulting 
from general expressions obtained in the previous sections. Finally, a brief conclusion is 
given in the last Sec. V. 



II. AVERAGE CURRENT FORMULA 



The study is started by considering a structure characterized by a sole lattice temperature, 
and a chemical potential which can be determined by the applied bias for each real contact. 
The current, in principle, is the expectation value of the current operator, which is generally 
written in the form 
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dv dr 

where \l/ CT (r, £) is the field operator with the spin-index a. Here, we neglect this index for 
simplicity and replace the sum by the factor i] s = 2 for the spin degeneracy. The transport 
is assumed to be along one direction so that we can expand the field operators in terms of 
transverse modes, {xo,n(ri)}, and of incident and outgoing waves, 

Xa,n(r±) 



[a a , n (t)e ika ' nZ + b a , n (t)e- ik ^ z ] (2) 
[tW^ + tWe^ 2 ], (3) 



where the index a is added to refer to the terminal the current is measured in; 
a ani a a,n), b^ n (b a ^ n ) are operators creating (annihilating) an electron in an incoming, outgo- 
ing state of the mode n, respectively; k a>n = ^j2m(E — E a ^ n )/h 2 , and v a)Tl is the velocity. 
Accordingly, the current operator (pQ) can be rewritten in the forrn^ 



T a ( z ,t) = 4X) / dEdE'e* E - E 'W h 

[al n {E)a a , n {E')-bl n {E)b a , n {E')]. 



(4) 



Note that Eq.(J3]) was obtained with the same assumptions as in the original standard 
scattering theory. It is now our key point of the analysis. Indeed, in the standard scattering 
theory one can set a linear relationship between the outgoing wave and the incoming wave 
amplitudes because of the phase coherence between them. However, in the presence of scat- 
tering processes such as the electron-phonon coupling, the phase of electrons is (randomly) 
changed when the particles jump into the scattering region. When they jump out, they 
have their own phase which is setup by the interacting processes they have experienced. 
Accordingly, the scattering region itself should be considered as a quasi-independent inter- 
mediate source (sink) providing (absorbing) electrons to (from) the collection (emission) 
terminal. So that in order to distinguish the coherent and incoherent particle fluxes through 
the structure we introduce a couple of operators d n and d^ n which characterize the processes 
of annihilating/ creating a particle in an interacting state in the scattering region. We thus 
propose the relation 
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where [s Q7jn z] still refers to the scattering matrix and [X ajT a] now refers to the dephasing 
matrix. The relationship among elements of these matrices is setup by the requirement that 
the operators ([5]) and ([6]) have to satisfy the anti-commutation algebra. That results in 
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Obviously, if operators d^ and d m are missing in Eq.(j7j) (X a = 0), the scattering matrix 
becomes unitary and the usual probability conservation, T n + R n = 1, is maintained. Under 



the current consideration, such unitary property is broken, but Eq.([7j) still leads to a more 
general conservation law as being seen later. We now denote 

Af^a, E, E') = 5 aS 5 ay 5 kl - ^aS,nk &a-y,nl 

n 
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C kl (a,E,E') = XXiJi^WJO. (8) 
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the current operator is then rewritten in the form 
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This equation obviously expresses all possible physical processes of particles such as the 
transmission from-lead-to-lead, from-lead-to-intermediate region-to-lead, and the relaxation 
itself. To get the expectation value of this operator it should be noticed that the average is 
now taken in quantum interacting states due to the presence of incoherent agents such as 
magnetic impurities, phonons, etc... However, we simply use 

(al n (E)a^ m (E')) = 8 a/3 S mn S(E - E')f a (E), (10) 
(al n (E)d m (E')) = 0, (11) 
(di(E)d m (E')) = 5 mn 5(E-E')f w (E), (12) 

where f a (E) is as usual the Fermi-Dirac function of electrons in the a-contact, and f w is 
referred to the occupation probability in interacting states in the device active region. Note 



that the appearance of S mn 5(E — E') in Eg. ([12]) does not mean only elastic processes being 
included, but it means /„, as the occupation probability of charges in excited states. The 
average of the second and third terms can be set to because they are equal in magnitude 
but with opposite sign, and therefore exactly cancel each other. After some simple algebra 
taking into account Eq.([7|) the current in the emitter writes 

7 «(*) = f E/ dE{[l-R etn (E)][f e (E)-f w (E)} 

n 

-T n (E)[f c (E)-f w (E)}}. (13) 

Here we have used the eigen-value representation of the matrices [s* ee (E)s ee {E)\ and 
[s* c (E)s ec (E)}, which correspond to the reflection (-R ej „) and transmission (T n ) probability of 
a transmission event through the whole structure in the conduction channel n. Note that we 
generally distinguish the reflection coefficients at each contact but, of course, a unique set of 
transmission coefficients {T n } governs the coherent transport through the whole structure.— 
A similar formula is also easily written down for the current measured at the collector. 
The conservation of current allows symmetrizing the current formula by evaluating the 
occupation factor f w in the scattering region, 

_ Ae/e + A cfc 

fw ~ A e + A c ' (i4j 



where A e ( c ) n = 1 — (R e ( c ),n + T n ) > 0. Substituting this result into Eq. ffTBl yields 

^ = lE / dET n\E) if e (E) - f c (E)\ , (15) 

n 

where Tj ot = T n + A Ae, " A A c '" regarded as the effective/total transmission coefficient. This 
result is really interesting because of its simple and well-known form of the Landauer formula. 
As expected, the total current is a sum of two contributions: the coherent (associated with 
T n ) and sequential (associated with the second term in T tot ) parts. Eq. ffl5|) leads to the 
conductance 



^ „ V A 6in + A c 

which is also coincident with that derived using other methods, such as the voltage-probe 
method,-^ and the NEGF method.^ 
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III. CURRENT NOISE POWER FORMULA 



We now go further to establish a formula for the current noise when both direct and 
indirect tunneling processes are taken into account. As usual we introduce the operator 
8I a (t) = I a (t) — (I a (t)), and define the correlation function S a p(t,t') between the currents 
in contact a and contact (3 as 

S a/3 (t,t') = (8I a (t)8I (t') + 8I p (t')8I a (t)). (17) 

The current noise power S a p(u) is then obtained as the Fourier transform of the time cor- 
relation function by noticing that S a p(t,t'), in fact, depends on t — t'. 

Using the formula of current operator, Eq.fjHJ), calculating the noise power leads to evalu- 
ate the expectation value of the products of two couples of creation and annihilation opera- 
tors in the interacting states. The products of operators of the same kind can be evaluated 
using the same algebra as that in the noninteracting case, i.e., 

8 aS 8 Pl 8 kn 8 lm 8{E 1 - E 4 )8(E 2 - E^f^E^l - f p {E 2 )\ 
+8 a p8 lS 8 kl 8 mn 8(E 1 - E 2 )8(E 3 - E A ) f a {E{) f\(E 3 ) , (18) 

and 

(4(£0^2)4,(£3K(£ 4 )> ^ 

8 kn 8 lm 8{E l - E i )8{E 2 - E 3 )f w {E x )[l - f w (E 2 )} 

+8 kl 8 mn 8(E 1 - E 2 )8(E 3 - E 4 ) f^Ex) f w (E 3 ) , (19) 

where f w is defined by Eq. ffT^l) . However, it is inappropriate to treat similarly the products of 
two couples of different operators since they manifest the correlation of charges transferring 
through different regions of the structure, the many-particle effects. Generally, we set 

(al Ctk (E 1 )di(E 2 )dl n (E 3 )as t n(E 4z )) = 8 aS 8 kn 8i m 

x8(E 1 - E 4 )8(E 2 - £ 3 )5 a (£i,£ 2 ), (20) 

and 

{d\ n {E 1 )a a ^i(E 2 )a\ m {E 3 )d n {E i )) = 8 a5 8 kn 8i m 

x8(E 1 - E 4 )8(E 2 - E 3 )E a (E x , E 2 ), (21) 
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where A a = S Q + S a implies the correlation of an occupied state in the contact a with an 
empty one in the scattering region, and vice versa [see Eq. pBl ]. In the result, we obtain a 
general formula for the noise power which may be expressed in four terms 

Sap(u) = + Sip + + S^p, (22) 

which are defined as 

e 2 f 

Sip = - V dE[A lS (a,E,E + hLu)A s ^(P,E + hLU,E)}' 
' nh<<-fj 

70 

x{f 7 (E)[l - f s (E + M + f s (E + M[l " ME)}, (23) 
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S «P = ^E dE[B 7 (a,E,E + fiw) 
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x 



B 7 (/3, £ + hiu, E)]'{E 7 {E, E + hio) + H 7 (£ + hu, E)}, (24) 



^ = ^E/ dE[B 7 (a, E,E + fko) 



x 



£ 7 (/3, E + fkj, E)]'{Ey(E + frw,E) + S 7 (£, E + hu)}, (25) 
e 1 r 

Sip = — dE[C(a,E,E + frw)C([3,E + hu;,E)]' x 
7th J 

{f w {E)[l - f w {E + hu) + f w (E + hu)[l - f w {E))}. (26) 

Here we use the notation [...]' for the trace of matrix in brackets. 

As in the standard scattering theory, although S a p is written at arbitrary frequency, 
equations (I22H26P are only applicable to the case of low frequency (for detailed discussions, 



see Ref. (4]). In this work, we are interested only in the zero-frequency noise for the two- 
terminal systems so that we can explicitly achieve a general expression, which reads 

S ee (0) = ^E fdE{T n [f e (l-f e )+f c (l-f c )}+T n (l-T n )(f e -f c ) 2 } 

n 

A 2 p 

+ X E / dE {KnW ~ fe) + /«,(! " fw)] + A e ,„(l - A e , n )A e } 
n 

+X E / dEA en T n [(f e - f c )(l - 2f e ) - (A e - Ac)]. (27) 

n 

It is clear that each term of Eq. (1271) can be associated with an explicit and simple physical 
meaning. Indeed, the first term, which coincides with the main results in Refs. pfl, arid 
the second one are the thermal and partition noise of the coherent and incoherent tunneling 
current, respectively, while the third term is interpreted as the correlation between the two 
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current components. Apart from this correlation, i.e., the correlation between the direct and 
indirect tunneling events, Eq. (l2"7j) also depends on A e(c ), the correlation among the charges 
transferred between the scattering region and the contacts, which is governed by the full 
statistics of the former region. On average, A e ( c ) can be well approximated by 

A e ( c ) = A° (c) = / e ( c )[l - f w ] + f w [l - fe(c)}, (28) 



and our result, Eq.( l271) . then reduces exactly to Eq.(18) in Ref.[17|, which was derived using 
the Keldysh's Green functions. 



IV. DISCUSSION 

In this section we shall discuss our main results, Eqs. fll5|) and (|27|) . derived in the pre- 
vious sections. It is certainly not necessary to do so for the current formula (fT5l) since it 
absolutely coincides with that obtained by different approaches.— 1 ^ So we shall discuss only 
the noise expression. Although Eq. (1271) has a simple form and has clear physical meaning 
for each term, it is difficult to go quantitatively in more detail because of some undetermined 
phenomenological parameters such as A e , A c , etc. However, in this section we shall qual- 
itatively discuss some important aspects of this formula. The discussion is for the typical 
double barrier resonant tunneling structure. After that, we illustrate this discussion using 
a concrete results obtained in the framework of a microscopic description. 

First of all, it is instructive to see a bit the microscopic picture of the electron interaction, 
e.g., the electron-phonon coupling, in the scattering region. The fact is that an inelastic 
interaction process may result in significant peaks corresponding to excited states in the 
density of states (DOS) picture (E + nhuo, where hu is the phonon energy, n is an integer 
number, and Eq is a bound state in the quantum well). The transport of electrons from the 
emitter to the collector can be then described in the following way. Even the level E is 
not available for the tunneling (lower the emitter conduction band bottom), electrons can 
jump into the well due to a presence of available excited states in the permitted energy 
interval. They then either relax to other states in the well before escaping to the collector, 
or coherently tunnel out from these excited states. Apparently, the height of excited peaks 
in the DOS determines the coherent tunneling probabilities {T n }, i.e., the higher excited 
peaks, the greater T n (see Fig. 1). The total current then, as pointed out in the previous 
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sections, can be seen as comprising of two components, sequential and coherent. 

Now, look at the third term in Eq. fl27j) . it expresses the correlation between the sequential 
and coherent components of current, and gives rise a negative contribution, 

8e 2 



J dEA e T n (fe ~ fc)(fe - U). (29) 



Here we have used Eq. (1281) . The term becomes important if the coherent tunneling prob- 
ability T n is large. Accordingly, if the inelastic interaction process is strong enough, it may 
cause a reduction of shot noise. While the current increases due to the opening of additional 
conduction channels, it is reasonable to expect a noise suppression. This point, comment 
(a), is contradictory to conclusions deduced from the voltage probe approach, 4 but it is 



commonly accepted (see Ref. |2l|, for instance) and will be illustrated later. 

Next, we discuss the case of very small T n , which corresponds to the tunneling induced 
by very low (insignificant) excited peaks of the DOS. In this limit the first and last terms 
in Eq. (1271) can be neglected and we keep only the second one, 

See(0) « -|- / dE {Knifed ~ fe) + fM ~ U)] + A e ,n(l " A e , n )A e } . (30) 

n 

The quantities A e and A c are then directly interpreted as the probability that charges tunnel 
through the first and second barriers, respectively, from the probability conservation law, 
A e ( c ) + R e ( c ) ~ 1. Consequently, they can be assumed as constants for simplicity. It should 
be kept in mind that in the interaction problems effects of many particles generally have to 
be properly considered. Here, the information on these effects is contained in the corrections 
to f w and to A e , A c . Eq. ( 1301) can be now estimated for zero temperature and high bias as 

A e + A c 



4e 2 A e A c 
S ee (0) oc — — - — - — 

and at the same time, 



- — c — + (1 - A e )-^ — -K 
A e + A c A c 



(31) 



2e A e A c 

I ^ ~rn r~ ■ ( 32 J 

h A e + A c 



The Fano factor therefore is equal to 



A e + A c v e! A c 

= 1 " A AeA ! + (1 - A e ) Ae + Ac M e . (33) 
A e + A c v e > A c K J 
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Here we have used the expression A e = A° + 5A e = A c /(A e + A c ) + 5A e . Apparently, if 
the correction 5A e is neglected, Eq. (1331) results in a Fano factor smaller than one (F < 1). 
However, if the many-particles effects are appropriately taken into account, the following 
condition can be satisfied: 

Under this condition, the Fano factor becomes greater than one, F > 1, i.e., the noise is 
enhanced. From Eqs. (1331) and (1341) . we can additionally deduce two following comments: 

(b) inelastic interaction processes may cause a weak enhancement of the Fano factor, and 

(c) it is easier to observe F > 1 in asymmetric structures wherein A e > A c . 

In order to illustrate comments (a), (b) and (c) we use a microscopic description for the 
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for the model of a GaAs 



process of inelastic electron-phonon interaction developed in Ref. 
quantum well embedded between two GaAs/AlGaAs layers. Fig. 1 plots the DOS and the 
retarded self-energy to show the formation of the excited states due to the interaction. The 
DOS exhibits a main resonant peak E Q and a significant phonon-emission peak E Q +ftuJo while 
the self-energy does the three excited states, one due to the phonon absorption E Q — hw , 
and two due to the phonon emission (the one with Eq + ftwo is strong and the other with 
E + 2hwo is weak) which are in accordance to the first and second phonon-assisted tunneling 
(PAT) regime. 

In Fig. 2 we plot two curves of the Fano factor which are extracted for the symmetric 
case, A e = A c , and the asymmetric one, A e > A c . Obviously, the figure shows a suppression 
of F in the first PAT regime and a weak enhancement in the second PAT one. The figure also 
shows distinguishably that the noise in the PAT regime is stronger suppressed and higher 
enhanced in the case of asymmetric compared to the symmetric one, that is consistent with 
the comment (c). 

Thus, we have discussed some aspects of the general noise expression derived in Sec. III. 
It was shown that the inelastic scattering process may induce either suppression or enhance- 
ment of the shot noise (at least in resonant tunneling structures) . This conclusion, which was 
then additionally discussed in the framework of a microscopic description, is not in agreement 
with those deduced by other approaches in literature. The matter is that in our considera- 
tion the role of many-particle effects is strongly emphasized, while in other phenomenological 
works it is not appropriately taken into account. It is important here to mention that the 
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super- Poissonian noise of F = 1.6, observed by Li et al. 23 in GaAs/AlcuGao.eAs-DBRTS, 
was ensured to come from the electron-phonon interaction process. This experimental data 
can be seen as a potential support to our discussion on the noise behavior. 

V. CONCLUSION 

In conclusion, we have developed the scattering theory in an appropriate manner such 
that it can be able to describe the electron transport through an intermediate region wherein 
interacting agents can strongly destroy the phase of incident electrons or make them relaxing 
their energy. Our key point here is based on the idea of considering the scattering region 
as an intermediate source/sink of particles, which maintains the phase correlation among 
the collision processes. Therefore, we can distinguish the coherent and incoherent particle 
fluxes which contribute to the total current independently. A derived formula of the current 
[Eq. (1151) ] recovers available results in the literature using different approaches. A new general 
formula for the current noise power [Eq. (1271) ] is also established, from which one can recover 
well-known results from the limits of purely coherent to purely sequential transport. A 
short discussion for the case of strongly incoherent tunneling in a double-barrier resonant 
tunneling structure is also presented and illustrated by a microscopic description. 
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FIGURE 1: Density of states and retarded self-energy describing the microscopic many- 
particle picture of the inelastic electron-phonon interaction. 

FIGURE 2: Fano factor curves extracted for the symmetric (Sym.) and asymmetric (Asym.) 
barriers of the double barrier resonant tunneling structures. 
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FIGURE 2 (V. NAM DO) 



16 



